Nuclear transparencies from photoninduced pion production 
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We present a relativistic and cross-section factorized framework for computing nuclear transparen- 
cies extracted from A(y, nN) reactions at intermediate energies. The proposed quantummechanical 
model adopts a relativistic extension to the multiple-scattering Glauber approximation to account 
for the final state interactions of the ejected nucleon and pion. The theoretical predictions are com- 
pared to the experimental 4 He(7,p7r~) data from Jefferson Lab. For those data, our results show 
no conclusive evidence for the onset of mechanisms related to color transparency. 
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The strong force exhibits a strong scale dependence. 
At low energies, hadrons are undoubtedly the adequate 
degrees of freedom. The properties of nuclei can be 
fairly well understood in a picture in which nucleons ex- 
change mesons. At high energies, that particular role 
of fermions interacting through force-carrying bosons, is 
played by quarks and gluons. The transition energy re- 
gion is a topic of current intensive research, for example 
at Jefferson Lab, where hadronic matter can be studied 
with intense beams of real and virtual photons possess- 
ing the proper wavelengths in the femtometer and sub- 
femtometer range. A commonly used observable to pin 
down the underlying dynamics of hadronic matter is the 
nuclear transparency to the transmission of hadrons. The 
nuclear transparency for a certain reaction process is de- 
fined as the ratio of the cross section per target nucleon 
to the one for a free nucleon. Accordingly, the trans- 
parency is a measure for the effect of the medium on the 
passage of energetic hadrons. It provides an excellent 
tool to search for deviations from predictions of mod- 
els based on traditional nuclear physics. One such phe- 
nomenon is color transparency (CT). Color transparency 
predicts the reduction of final state interactions (FSI) of 
hadrons propagating through nuclear matter in processes 
at high momentum transfer. Experiments have been car- 
ried out to measure nuclear transparencies in search of 
CT in A(p, 2p) 0, S, 1 1 and Me^v) U S S & i M 
reactions, p-meson production |lllll2i| and diffractive dis- 
sociation of pions into di-jets [13j |. Intuitively one expects 
CT to reveal itself more rapidly in reactions involving 
mesons. Indeed, it appears more probable to produce a 
two-quark configuration with a small transverse size and, 
as a consequence, reduced FSI. 

Recently, the nuclear transparency for the pion photo- 
production process 771 — > 7r~_p in 4 He has been measured 
in Hall A at Jefferson Lab [14]]. In the experiment, both 
the outgoing proton and the pion are detected for pho- 



ton beam energies in the range 1.6-4.5 GeV and pion 
center-of-mass scattering angles of 70° and 90°. With 
two hadrons in the final state, the pion photoproduction 
reaction on helium offers good opportunities to search 
for medium effects that go beyond traditional nuclear 
physics. Not only is He a dense system, its small ra- 
dius optimizes the ratio of the hadron formation length 
to the size of the system. Along the same line, a pion 
electroproduction experiment at Jefferson Lab, measur- 
ing pion transparencies in H, D, 12 C, 64 Cu and 197 Au [l5| 
has been completed and data analysis is currently under 
way. It speaks for itself, though, that the availability 
of a model for computing transparencies extracted from 
A(7, 7riV) and A(e, e'nN), is essential for interpreting all 
these measurements. As a matter of fact, it is a real 
challenge to compute the effect of the medium on the si- 
multaneous emission of a nucleon and a pion within the 
context of traditional nuclear-physics models. 

In Ref. [l4[ the results for 771 — » ir~p on 4 He are 
compared to the predictions of a semi-classical model by 
Gao, Holt and Pandharipande [lj|. The results of the 
calculations including CT effects were found to be more 
consistent with the measurements. An alternative semi- 
classical approximation to compute pionic transparencies 
is presented in (l7| . A quantummechanical model for me- 
son photo pro duction on 2 H was recently proposed by J.- 
M. Laget [18| . In this Letter, we present a quantumme- 
chanical model for computing the nuclear transparencies 
in 771 — ► Ti~p on finite nuclei. The model is based on 
a relativized version of Glauber theory and is essentially 
parameter-free. To our knowledge the presented frame- 
work is the first of its kind. 



In describing the ^(7, Ntt)A — 1 reaction we use 
the following lab four- momenta: Q tl (q,q) for the 
photon, P A (E A ,p A 
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O) for the target nucleus, 
P A -i(Ea-i,1>A-i) for the residual nucleus, Pfc — 
(En,Pn) and P% = (E n ,p n ) for the ejected nucleon 
and pion. The missing momentum p m is defined as 
p rn = —pA-i = pn + Pit — <f- Following the conventions 
of Ref. 19], the fivefold differential cross section reads in 
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the lab frame 
da 



dE n dQ,.„dn N 4:(2n) 5 qE 



with the recoil factor given by: 
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and M^' 1 ^^ the invariant matrix element: 



= (P£,P»m s ,PZ_ 1 J R M R \6\Q»,PZ0 + ) , (3) 



where m s is the spin of the ejected nucleon N and JrAIr 
the quantum numbers of the residual nucleus. We restrict 
ourselves to processes with an even-even target nucleus 
A. The operator O describes the pion photoproduction 
process and we assume it to be free from medium effects. 
This is a common assumption in nuclear and hadronic 
physics and is usually referred to as the impulse approx- 
imation (IA). 

For the target and residual nucleus we use relativis- 
tic wave functions as they are obtained in the Hartree 
approximation to the crw-model with the Wl parameter- 
ization [20] . When studying the transparency, it is conve- 
nient to factorize the invariant matrix element M^2' N ^ 
into a part containing the elementary pion photoproduc- 
tion process and a part with the typical medium mecha- 
nisms in the process under study. It is clear that the at- 
tenuation on the ejected proton and pion induced by FSI 
mechanisms belongs to the last category and determines 
the nuclear transparency for the process under study. 
Even in the relativistic plane-wave limit for the ejected 
nucleon and pion wave function, factorization of the cross 
section is not reached through the presence of negative- 
energy contributions. Neglecting these, the computation 
leads to an expression for the cross section in the rela- 
tivistic plane wave impulse approximation (RPWIA) 
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with p a {p m ) = J2 ms , m \ u (Prn,m s )(f) a (p m )\ 2 the momen- 
tum distribution which is obtained by contracting the 
bound-state wave function <p a with the Dirac spinor u. 
The a(n,K,m) denotes the quantum numbers of the 
bound nucleon on which the photon is absorbed. Further, 
denotes the cross section for 7 + N — ► it + N', and 

and t = (Q>*- P%f are the Mandelstam 
variables . 

In this work, we concentrate on A^, Ntt)A — 1 pro- 
cesses for which the wavelengths of the ejected nucle- 
ons and pions are typically smaller than their interac- 
tion ranges with the nucleons in the rest nucleus. Those 
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conditions make it possible to describe the FSI mecha- 
nisms with the aid of a Glauber model. A relativistic 
extension of the Glauber model, dubbed the Relativistic 
Multiple-Scattering Glauber Approximation (RMSGA), 
was introduced in Ref. [21]. In the RMSGA, the wave 
function for the ejected nucleon and pion is a convolu- 
tion of a relativistic plane wave and an eikonal Glauber 
phase operator S-psiir) which accounts for all FSI mech- 
anisms. Through the operation of <Sfsi(^) every residual 
nucleon in the forward path of the outgoing pion and nu- 
cleon adds an extra phase to their wave function. The 
RMSGA framework has proved succesful in describing 
cross sections and other observables in exclusive A(e, e'p) 
[21], [22j and A(p, 2p) [23| reactions. The numerically chal- 
lenging component in RMSGA is that <Sfsi(^) involves a 
multiple integral which tracks the effect of all collisions 
of an energetic nucleon and pion with the remaining nu- 
cleons in the target nucleus. Realistic nuclear wave func- 
tions arc also used in the models of Refs. [lflQjl- Con- 
trary to the RMSGA model, however, the transparencies 
are computed at the squared amplitude level adopting a 
semi-classical picture for the FSI mechanisms. 

In the numerical calculations within the context of the 
RMSGA, the following phase is added to the product 
wave function for the ejected nucleon and pion: 



SfsiW = n [l - IW(?- ti)9{zj - z) 

3=2 

x [l-T nN (b' -btyiz'j-z') 



(5) 



where fj(bj, Zj) are the coordinates of the residual nucle- 
ons and f(b, z) specifies the interaction point with the 
photon. In Eq. ([5]), the z and z' axis lie along the 
path of the ejected nucleon and pion respectively. The 
b and b' are perpendicular to these paths. Reflecting 
the diffractive nature of the nucleon-nucleon (N'N) and 
pion-nucleon (irN) collisions at intermediate energies, the 
profile functions Tn>n and T^n hi Eq. ([5]) are parame- 
terized as 
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(with, i = ir or N') 
(6) 



Here, the parameters a\^ (total cross section), [3im (slope 
parameter) and CiN (real to imaginary part ratio of the 
amplitude) depend on the momentum of the outgoing 
particle i. In our calculations those parameters are 
obtained by interpolating data from the databases for 
N'N — ► N'N from the Particle Data Group [H and 
irN — ► irN from the analysis of Refs. [25l l2q. 

Now we derive an expression for the fivefold 
A(7, Nir)A — 1 cross sections when implementing FSI ef- 
fects. To this end, we define the distorted momentum 
distribution: 
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Here, tf>%(p) = (2 ^ 3/2 / dre ^^(^^(f) is the dis- 
torted momentum-space wave function, which is the 
Fourier transform of the bound nucleon wave function 
and the total Glauber phase. In the absence of FSI, 
the pfiMSGA^Pm) of Eq. ([7|) reduces to the p a {p m ) in 
Eq. ^ when negative-energy components are neglected. 
Based on this analogy we obtain the cross section in the 
RMSGA approach by replacing p a {p m ) by P RM sga(p™) 
in Eq. ©. 

In our calculations, color transparency effects are im- 
plemented in the standard fashion by replacing the to- 
tal cross sections cr-^ in Eq. © with effective ones [27| 
which account for some reduced interaction over a typical 
length scale lh corresponding with the hadron formation 
length (i — ir or N') 
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Here n is the number of elementary fields (2 for the 
pion, 3 for the nucleon), k t — 0.350 GeV/c is the av- 
erage transverse momentum of a quark inside a hadron, 
Z is the distance the object has travelled since its cre- 
ation and lh — 2p/ AM 2 is the hadronic expansion length, 
with p the momentum of the final hadron and AM 2 the 
mass squared difference between the intermediate pre- 
hadron and the final hadron state. We adopted the values 
AM 2 = 1 (GeV) 2 for the proton and AM 2 = 0.7 (GeV) 2 
for the pion. 
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FIG. 1: The nuclear transparency extracted from 4 He(7,p7r~) 
versus the squared momentum transfer t | at — 70°. 

The solid (dashed) curve is the result of the RMSGA calcu- 
lations without (with) color transparency. The semi-classical 
model [l6|] results are presented by the shaded areas: the 
hatched (dotted) area is a calculation without (with) CT. 
Data from 0. 

In Figs. Q] and [2 we present the results of trans- 
parency calculations for He together with the exper- 
imental data and the predictions of the semi-classical 
model of Ref . [l~6| . In comparing transparency measure- 
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FIG. 2: As in Figure [TJ but now for 6£ m . = 90°. 



ments with theory, accurate modeling of the experimen- 
tal cuts is required. We adopt the following definition for 
the transparency 



T = 



RMSGA 



dE^.dn^.dnN. I 



(9) 
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where i denotes an event within the ranges set by the 
detector acceptances and applied cuts. Further, ^ Q ex- 
tends over all occupied single-particle states in the target 
nucleus. All cross sections are computed in the lab frame. 
Further, Y(q) is the yield of the reconstructed experimen- 
tal photon beam spectrum for a certain photon energy 
[i~4j . We assume that the elementary 7 + n — > n~ + p 
cross section d ^ in Eq. (fj| remains constant over the 
kinematical ranges which define a particular data point. 



With this assumption the cross section 
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cancels out 



of the ratio ([9]). In order to reach convergence in the 
phase-space averaging in Eq. @ we generated about 
one thousand theoretical events within the kinematical 
ranges of the experimental acceptances. This was done 
for all data points, eight in total, and corresponding kine- 
matical ranges, of the Jefferson Lab experiment. Detailed 
kinematics for these data points can be found in Ref. [3] • 
The computed RMSGA nuclear transparencies are sys- 
tematically about 10% larger than the ones obtained in 
the semi-classical model. As can be seen in Fig. [TJ our 
model predicts a rise in the transparency for | t \ val- 
ues below 1.2 GeV 2 . This rise is due to the minimum 
in the total proton-nucleon cross section in Eq. ^ for 
the proton momenta associated with these momentum 
transfers. The RMSGA results overestimate the mea- 
sured transparencies at small | t |, but do reasonably 
well for the higher values of | t \. Inclusion of CT effects 
tends to increase the predicted transparency at a rate 
which depends on a hard-scale parameter. Here, that 
role is played by the momentum-transfer | t |. Thus, in- 
clusion of CT mechanisms results in an increase of the 



4 



nuclear transparency which grows with the momentum 
transfer | t |. The magnitude of the increase depends on 
the choice of the parameters in Eq. (JU). For the moment, 
there are no experimental constraints on their magnitude. 
As can be appreciated from Figs. Q] and [2l the RMSGA 
calculations predict comparable CT effects as the semi- 
classical calculations. We have to stress though that the 
calculations with CT are normalized to the calculations 
without CT for the data point with the lowest \t\ in the 
semi-classical model. We did not perform this normal- 
ization for our calculations. Our results without color 
transparency are in better agreement with the experi- 
mental results than those with CT effects included. This 
is in disagreement with the semi-classical model whose 
results with CT effects are in better agreement with the 
experimental data. We also have to point out that, al- 
though the calculations with CT effects overestimate the 
experimental results for all data points, the slope of this 
curve shows better agreement with the slope of the data 
than the slope of the curves without CT effects. 

To provide an idea of the A-dependence of the nu- 
clear transparency extracted from A(7,p7r~), we plotted 
in Fig. [3] the calculations for one data point for several 
nuclei with the same kinematic cuts as before. However, 
due to these cuts, nucleon knockout from the innermost 
shells in the heavier nuclei was not always possible for the 
generated events in the calculations. The transparency 
would even be lower for these heavier nuclei if no cuts 
would be applied. 



3 0.6 



0.5 



0.3- 



0.2 



FIG. 3: A-dependence of the transparency. Calculations were 
made for 4 He, 12 C, 16 0, 56 Fe and 208 Pb at |t| = 3.5 GeV 2 for 
#c.m. — 90° without color transparency 



In summary, we have developed a quantum mechani- 
cal model based on a relativistic extension to multiple- 
scattering Glauber theory to calculate nuclear trans- 
parencies extracted from ^(7, iV7r) processes. The model 
can be applied to any even-even target nucleus with a 
mass number A > 4. The nuclear transparency is the 
result of the attenuating effect of the medium on the 
ejected proton and pion, and is computed by means of 
a Glauber phase operator. The numerical computation 
of the latter, requires knowledge about irN — > irN and 
N'N — > N'N cross sections, as well as a set of rela- 
tivistic mean-field wave functions for the residual nu- 
cleus. In contrast to alternative models, which adopt 
a semi-classical approach, we treat FSI mechanisms at 
the amplitude level in a quantum mechanical and rela- 
tivistic manner. Comparison with experimental results 
for helium shows no evidence of color transparency in 
our model. Further progress will very much depend on 
the availability of new data. The model presented here 
can be readily extended to electroproduction processes 
for comparison with the forthcoming Jefferson Lab data. 
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